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SUMMARY - We present here the basic issues of our simulations of Planck 
observations, focusing on the study of the stripes generated by l//-type noise due 
to amplifiers noise temperature variation of LFI radiometers generated by gain 
fiuctuations. Our simulations include a realistic estimate of 1// knee frequency, 
based on a recent analytical study of the properties of the Planck LFI radiometers 
and of their systematic effects, as well as realistic choices of load, amplifier noise 
and paylod environment temperatures. By comparing simulated observed maps 
of CMB anisotropics obtained by including or not the instrumental noises we are 
able to quantify the magnitude of the striping effect. The "standard" destriping 
method presented in the report on the Phase A study has been adapted to our 
purposes. We quantify the efficiency of this destriping technique for different 
scanning strategies, for typical off-axis LFI beams. Our results are compared 
with those of other previous analyses. 

1 Introduction 

The Planck Surveyor is an European Space Agency (ESA) satellite mission to map spatial 
anisotropy in the Cosmic Microwave Background (CMB) over a wide range of frequencies with 
an unprecedented combination of sensitivity, angular resolution, and sky coverage (Bersanelli 
et al. 1996). The data gathered by this mission will revolutionize modern cosmology by 
shedding light on fundamental cosmological questions such as the age and present expansion 
rate of the universe, the average density of the universe, the amount and the kind of dark 
matter, and other questions. As with any CMB experiment, achieving the desired perfor- 
mance requires careful attention to the control of systematic effects. l//-type noise in the 



radiometer output is one of the most critical systematic effect pertaining to the Low Fre- 
quency Instrument (LFI) radiometers because it may lead to striping in the final sky maps 
and increase the noise level. In general, a value of the 1// knee frequency, fk, significantly 
greater than the spacecraft rotation frequency, /«, will lead to some degradation in sensitivity 
(Janssen et al. 1996). In this paper we examine the impact of the l//-type noise by adopting 
a realistic estimate of 1// knee frequency as a function of the load, amplifier noise and pay- 
lod environment temperatures, of the radiometer bandwidth and of the level of fluctuations 
of each stage of the radiometer ampliflers, based on recent analytical studies of systematic 
effects in Planck LFI radiometers (Seiffert et al. 1997). 

In section ^ we summarize the basic concepts relevant for understanding the basic prop- 
erties of instrumental noises. 

In section ^ we present the mathematical formalism of our numerical code for the sim- 
ulation of the Planck observations, including instrumental noises, and the data stream gen- 
eration; we typically refer here to Planck like scanning strategies, but our code is versatile 
enough to allow the study of other observational schemes. 

In section ^ we describe how we have converted our simulated data streams in simulated 
observed maps. 

In section || we discuss in detail the mathematical formalism of the proposed destriping 
technique, including some considerations about the numerical efficiency. 

The "standard" estimators that quantify the magnitude of the striping effect and the 
efficiency of destriping techniques are presented in section g. 

Some preliminary results are presented in section 1^. 

Finally, in section 8 we draw out the main conclusions of our analysis, compare our results 
with those of previous works and draw out a brief guide-line for a future work. We discuss 
there the main implications of our study, focusing on their impact for the optimization of the 
Planck observational strategy. 

2 Sources of instrumental noises 

Planck LFI radiometers are modified Blum correlation receivers (Blum 1959, Colvin 1961). 
The modification is that the temperature of the reference load is quite different from the sky 
temperature (Bersanelli et al. 1995). To compensate, differing DC gains are applied after 
the two detector diodes. Adjusting the ratio of DC gains, r, allows one to null the output 
signal, minimize sensitivity to RF gain fiuctuations, and achieve the lowest white noise in 
the output. Although it may not be immediately apparent, the fact that the reference load 
is not at the same temperature as the sky does not increase the white noise level compared 
to a standard correlation receiver. 

The ideal sensitivity of our radiometer for a single observation with an integration time 
r is 



AT- \/2 {Tn + Tr) , , 

^J white = 7^= , Uj 

where P is the effective bandwidth, T^ is the noise temperature of the signal entering one of 
the two horns and Tn is the amplifier noise temperature. 

In order to null the average output signal of the radiometer, the DC gains ratio after 
the two detector diodes, r, must be adjusted to the proper value. In the simple case that 
the gains of the two amplifiers entering the two horns and their noise temperatures can be 
considered equal, r must be setted to the value: 



r = TfT—TfT , (2) 

where Ty is the reference load temperature entering the other horn. The above temperatures 
are antenna temperatures; T^ is due properly to the sum of the sky temperature (essentially 
the CMB monopole antenna temperature, related to the CMB thermodynamic temperature 
To ~ 2.726K, plus "minor" contributions from CMB dipole and anisotropies, galactic and 
extragalactic foregrounds, bright sources, Zodiacal light ...) and of the "environment" tem- 
perature (of about 1 K) due to the satellite emission. 

There are several potential concerns for the current radiometer scheme. Amplifiers noise 
temperature variations, that derive from gain fluctuations, could be confused with the sky 
signal variations, that we are interested in measuring, by introducing a change in the observed 
signal, ATequiv, which can mimic a true sky fluctuation. The amplifiers noise temperature 
variations have the characteristics of 1// noise and this leads to l//-type noise in the ra- 
diometer output. 

We have recently estimated the 1// knee frequency of our radiometer (Burigana et al. 
1997a, Mandolesi et al. 1997a, Seiffert et al. 1997) on the simplifying assumption that there 
is no 1// contribution from the detector diodes and neglecting the effects of the phase shifting 
that is designed to control this contribution; the diodes are assumed to be perfect square law 
detectors and we have assumed that the bandpass of the signal is the same in both legs of 
the radiometer. Here we briefly summarize the main arguments relevant for the expected 
magnitude of gain and noise temperature fluctuations and the main concepts relevant for the 
estimate of the 1// knee frequency. 

The contribution to the 1// noise directly due to amplifier gain fluctuations results to 
be zero at first order. Under quite reasonable assumptions, the noise contributions due to 
reference load fluctuations and fluctuations in the ratio of DC gains are much less than that 
due to the amplifier noise temperature fluctuations term, which then results to be the domi- 
nant source of 1// noise in the radiometer output. Imperfect isolation does not significantly 
modify this conclusion. Further, the sensitivity of our radiometer to differences between the 
gains and noise temperatures of the two amplifiers is not critical. As a consequence, all these 
complications cannot significantly change the knee frequency respect to the results we draw 
out here below. 

Cryogenic HEMT amplifiers have noise temperature fiuctuations with a l//-type spec- 
trum because induced by the l//-type gain fiuctuations of amplifiers (Wollack 1995, Jarosik 
1996, Seiffert et al. 1996). The magnitude of noise temperature fluctuations can be computed 
from the following argument. Assuming that each stage of the amplifier has the same level of 
fiuctuation, we can conclude that the transconductance of an individual HEMT device also 
fiuctuates according to: 

Ag^ _ 1 AG 



gm 2./Ws G 

where Ng is the number of stages of the amplifier, typically ~ 5. An optimal low noise 
amplifier design will have equal noise contributions from the gate and drain of the HEMT, 
which means the changes in g^ will lead to changes in T„ (Pospieszalsky 1989). This can be 
expressed as 

AT„ Agm ,.. 



We can write the 1// spectrum of the gain fluctuations as: 

— - — (5) 

Putting this together we get: 

We can therefore write the amphfier noise temperature fluctuations as 

AT„, A 



Tn ^/f' 



(7) 



with A = C/{2y/Ns); a normalization of A ~ 1.8 x 10^^ (relying on the references above) is 
appropriate for the Planck radiometers at 30 and 45 GHz. Throughout, we will use units of 
K/ vHz for AT so that we will not need to refer to the sampling frequency of the radiometer. 
In these units then, AT/T has units of Hz~^ " and A is dimensionless. We also note that the 
value of A will generally depend on the physical temperature of the amplifier. The values for 
A given here should be regarded as estimates rather than precise values. For the radiometers 
at higher frequencies, it will be necessary to use HEMT devices with a smaller gate width 
to achieve the lowest amplifier noise figure. We expect that the gate widths will be roughly 
1/2 that of the devices used for the lower frequency radiometers and this will lead to gm 
fluctuations that are roughly a factor of v2 higher (Gaier 1997, Weinreb 1997). We will 
therefore adopt a normalization of A = 2.5 x 10~^ for the 70 and 100 GHz radiometers. 

Starting from the expression of the average output of the differential radiometer, one 
derives the change in the output signal for the above small change in the noise temperature 
of one of the amplifier; and by multipling it by a factor \/2 because both amplifiers (which 
have uncorrelated noise) can contribute to this effect, we have the change in the observed 
signal, ATequiv, given by: 



^-'cauiv — V ^^J-n 



'- cquiv 



(8) 



We define the "knee" frequency as the post-detection frequency, fk at which the 1// 
contribution and the ideal white noise contribution are equal, i.e. ATgquiv = AT^hitc. For 
the computation of the knee frequency we use an integration time r = 1/(2 A/) and A/ = 1 
Hz, according to the choosen units for A and AT„. The knee frequency is then given by: 



This expression shows as the knee frequency depends upon several factors including the 
radiometer bandwidth, reference load temperature, and the intrinsic level of fluctuation in 
the HEMT devices; values of ~ few x 0.1 Hz, according to the frequency, should be reached 
with only passive cooling of the radiometer (to about 50 K), whereas active cooling (to about 
20K or less) can further reduces the knee frequency. As examples, assuming a 20% bandwidth 
for our frequency channels and an antenna temperature T^ = 3K, the knee frequency is 0.046 
Hz and 0.11 Hz respectively at 30 GHz (assuming Ty = 20K and T„ = lOK) and at 100 GHz 
(assuming Ty = 20K and r„ = 40K). 



For a more realistic evaluation we have to accurately repeat the analytical considerations 
by Seiffert et al. (1997) to distinguish between the amplifier noise temperature, r„, and the 
system temperature, Tgyg. In reality, it would be more appropriate to insert Tgys in eq. (1) 
and to carefully consider when the two temperatures enter in the determination of the knee 
frequency. On the other hand, the difference between these two temperature is estimated to 
be of few K while the instrumental situation is still partially unclear; then a careful distinction 
between them is presently not necessary from practical point of view, although interesting. 
We will address this point in a future work. 

The knee frequency must be compared to the spin frequency fs] for the Planck observa- 
tional strategy proposed for the Phase A study /s = 1 r.p.m., i.e. 0.017 Hz. 

As a comparison, for a total power radiometer, laboratory measurements have found knee 
frequencies between 10 and 100 Hz; the modified correlation radiometer scheme reduces the 
knee frequency by more than two order magnitudes. 

3 Simulation of the mission 

We have written a code that simulates the basic properties of Planck observations in order to 
study the striping effect due to l//-type noise on the measured sky temperatures. For sim- 
plicity our sky map includes only CMB fluctuations, generated in equi-cylindrical pixelisation 
(ECP) by using the method of Muciaccia et al. (1997); then we have projected it in COBE- 
cube pixelisation in order to have quasi equal area pixels. For our tests we have considered 
CMB fluctuations in a typical CDM scenario {Qh = 0.05, with COBE/DMi? normalization). 
We have used here an input map Min at a resolution of 19.4' [i.e. at COBE-cube resolution 9] 
that we have derived from an ECP map with 1024 grid points along each parallel computed 
by using the multipoles up I = 512. The typical dimension of a pixel of an input map at 
resolution 9 is comparable with the beam FWHM of 30' of the channels at 30GHz, that we 
have considered for the present tests. Our input map is shown in Figure 1. 

3.1 The simulation of the sky observation 

The instrument design for Planck mission calls for multi-frequency focal plane arrays placed at 
the focus of off-axis optical systems, in order to achieve proper angular resolution, sensitivity, 
and spectral coverage. As a consequence, not all the feedhorns can be located very close 
to the centre of the focal plane; the study of the implications of the related beam optical 
distortions on temperature measurements has been presented in other works (Burigana et al. 
1997b, c). For the present purpose we do not convolve the sky map with the beam response, 
but we simply read the map temperature corresponding to the pixel identified by the central 
direction of a given beam during the sky scanning. 

Figure 2 shows the schematic representation of the observational geometry. 

Let i be the angle between a unit vector s, along the satellite spin axis (outward the 
Sun direction) , and the normal to the ecliptic plane, and p the unit vector of the direction 
of the optical axis of the telescope, at an angle a from the spin axis (i = 90° and a = 70° 
for the Phase A study (Bersanelli et al. 1996). We choose two coordinates x and y on the 
plane tangent to the celestial sphere in the telescope optical axis direction, with unit vector 
u and V respectively; we choose the x axis according to the condition that the unit vector 
u points always toward the satellite spin axis; indeed, for standard Planck observational 
strategy, this condition is preserved as the telescope scans different sky regions. With this 
choice of reference frame, we have that v = pA s/\pA s\ and u = vAp/\vAp\ (here A indicates 
the vector product). In general the coordinates (xo,yo) of the beam centre will be identified 
by two angles; we use here the colatitude 6b and the longitude (f>B in the u, v, p reference 



frame (see the Appendix A for details on geometrical transformations). For the present test 
we assume a typical off-axis location of the considered beam: 6b = 2.8°, cpB = 45°. We 
note that our choice of 9b is representative of typical LFI beam position for a telescope with 
a primary mirror of 1.5m aperture in the new Planck optical configuration (see Mandolesi 
et al. 1997b). Our choice of (pB corresponds to a case of intermediate efficiency respect 
to the destriping technique; the cases (pB = 0° (beam located along the u direction in the 
TICRA U-V plane) and (pB = 90° (beam located along the v direction in the U-V plane) 
are equivalent respectively to an on-axis beam and to a beam that suitably distributes the 
crossings between different circles in two regions, close to the ecliptic poles, of maximum 
size. The off- axis choice for 9b, when cpB significantly differs from 0, also assures that, even 
in the particular case of an angle of 90° between the spin axis and the telescope direction, 
the scanning circles do not cross always exactly the ecliptic poles but two somewhat larger 
regions around them also for scanning circles corresponding to significantly different spin axis 
positions. 

We remember that in the Planck scanning strategy of the Phase A study the sampling 
time of each receiver is choosen in order to have three samplings when the telescope axis 
describes in the celestial sphere an angle with a length equal to the beam FWHM (~ three 
samplings per beam). For any value of the angle a, this condition determines the number of 
samplings, Up, per scan circle. 

We will study in detail the effect introduced by the telescope motion in a future work. For 
the present purpose, we simply assign the proper beam directions to each sampling, extract 
the corresponding pixel and read its temperature in the simulated input map; given the 
assumed FWHM of 30' and the pixel dimension 19.4' we have typically 2-3 different pixels 
explored in an integration time corresponding to 3 samplings. 

We remember that, in order to be able to reduce the white noise in a simple way, we 
want to "close" the scanning circle, i.e. we need to modify just a little the integration time 
by requiring that the telescope points always at the same set of directions when it repeats 
the (120 for the Phase A study) number of cycles with the same spin axis direction. Finally 
we have applied a further reassessment of the integration time to make Up multiple of 12, in 
order to be able of applying in the future the destriping technique not only with one but also 
with 2, 3 or 4 level constants for circle (see sections |5| and \l^ ). 

Our code let us free to implement arbitrary scanning strategies (see also Appendix A); 
we consider here cases with the spin axis always on the ecliptic plane and with an angle a 
between the spin axis and the telescope direction of 80°, 85°, 90°; lower values of a, like that 
of 70° assumed for the Phase A study, will require wide oscillations (15° — 20°) of the spin axis 
on the ecliptic plane (with relevant problems for the thermal stability) in order to observe 
the regions close to ecliptic poles, that are indeed very informative, being not significantly 
contaminated by the galactic emission. We have considered also a case with a = 90° and 
with ten 10° sinusoidal oscillations of the spin axis on the ecliptic plane. For all the cases 
we have considered a 360 days mission in order to compare results obtained with the same 
mission duration. 

3.2 The generation of instrumental noise series 

We generate white noise and 1/f noise using a random number generator code (see Press et al. 
1992). For the white noise this is very simple: we have simply to rescale a gaussian random 
noise distribution, by taking into account the radiometer rms white noise [see eq. |l|. For 
including the l//-type noise we have adapted in FORTRAN an original IDL code provided us 
by M. Seiffert, based on the power spectrum expansion in the Fourier space of the noise com- 
ponents. We generate together white and 1/f noise. Firstly a random gaussian distribution 



of (white) noises is generated; then we calculate its power spectrum using an FFT code (see 
Press et al. 1992); the amplitude of this power spectrum is then multiplied by (1 + fk/ fY 
for including the 1// contribution. Finally we have the time series with both the noises by 
computing the (inverse) FFT of the power spectrum. When we take into account the real 
properties of the 30 GHz receivers, eq. (9) (Seiffert et al. 1997) gives the 1// knee frequency 
of our radiometers. We adopt as reference for the present simulations fk = 0.05Hz a value 
adequate to a cooling efficiency that allows to keep a load temperature Ty ~ 20K (we will use 
here r„ = 9K, /3 = 6 GHz). Our noise series consist of about 2 x 10^ evaluations and a single 
series covers the integrations for 8 different contiguous spin axis directions (16 hours). We use 
a single series for white noise and 1// noise together, but we generate also a noise series that 
includes white noise only; this increases the computation time of only few percent (see above) 
but it is useful for comparison and for the quantification of the stripes magnitude and of the 
destriping efficiency (see sections 6-^ 10). In our code the generation of noise series is coupled 
with the sky observation; more precisely, because of the noises dependence on the antenna 
temperature Tx, we include the exact local sky temperature in the noise magnitudes. This 
kind of flexibility may be useful also in the future for simulations that will include Galaxy 
emission and thermal drifts too. Then, the only simplification is in the estimate of the knee 
frequency for the 1// noise series, which is assumed to be constant (i.e. with constant Tx, 
not allowing for spatial or time variations); the adopted value is consistent with the CMB 
monopole antenna temperature and a typical environment temperature of 1 K. 

3.3 The data stream generation and recording 

We have recorded our data streams in 4 matrices; any row of these matrices refers to the 
data obtained from a given spin axis direction; the number of columns is equal to the number 
of samplings (n^ ~ 2100, depending on a); the number of rows is equal to the number, n^, 
of different spin axis directions (4320 for a shift of 5' in the spin axis direction). We have 
recorded the following data: 

• the matrix N which contains the pixel numbers - 4 bytes per pixel - (in COBE-cube 
pixelisation) corresponding to the different integrations; 

• the matrix T which contains the "global" temperatures observed by the receiver in the 
above sky directions, directly averaged over the number of cycles per scan circle (120 
in the Phase A study scheme) in order to avoid the use of a large useless amount of 
memory space. T includes: the "input" sky temperature fluctuation as observed in the 
adopted geometrical scheme, which obviously can be computed a single time for any 
given spin axis direction; the white noise and the 1/f noise averaged over the number 
of cycles for each spin axis direction; 

• the matrix W, generated in the same way as the matrix T, but containing the temper- 
atures that will be observed in presence of the white noise only (see section 3.2); 

• the matrix G contains the temperatures that will be observed in absence of instrumental 
noise and will be useful to check the goodness of the geometrical part of our flight 
simulation code. 

Of course the matrices W and G do not have the corresponding ones in a real observation. 

In principle, if we simulate the mission for a time long enough that the spin axis return 
on the same direction after a certain period (360 days, as a reference allowing for the case 
a = 90°) we can average the data of the second period with the corresponding rows of 
the flrst period, and take memory that their will be affect by a (statistical) error v2 times 



smaller than that corresponding to rows that are observed for a single period. For a more 
realistic simulation, we need to record a further matrix E with the statistical sensitivities 
corresponding to the pixels of matrix N, which takes into account this fact as well as a 
possible degradation in sensitivity for some elements of data streams due for example to 
cosmic rays, spurious effects ... . We neglect these kinds of complications for the present 
analysis, which is equivalent to say that all the elements of E assume a constant value. 
Nevertheless, in section ^ we draw out our formulas by including this possible effect, for sake 
of generality. 

For implementing the destriping techniques (see section 5) we need to recognize when the 
pointing directions for different spin axis directions are substantially identical. The sky point- 
ing direction is stored in this scheme only through the corresponding pixel number. Then, 
from a statistical point of view, two pointing directions are considered identical provided that 
their distance is smaller than the pixel size. The resulting number of pixels in common is 
then related to the assumed pixel size. We have then recorded: 

• additional matrices Nh (H=l,2,...,n^) contains the pixel number corresponding to the 
different integrations for a certain number, n/j, of resolutions higher than that used for 
the input /output maps (for example at resolution 10 or 11 for input map at resolution 
9). By exploiting these matrices we will able to test the adopted destriping technique 
under more stringent conditions on the average distance between their pointing direc- 
tions in the research of the pixels in common. 

4 From data streams to observed simulated maps 

Given the above simulated data streams, it is quite simple to obtain the following simulated 
observed maps (see also section [5^ ) , that can be easily compared one to each other. 

We compute the sensitivity map, Ms, with which any map pixel is observed, by recogniz- 
ing how many times a given pixel is observed from the analysis of the whole matrix N (and 
E if this is the case) . 

From the matrices N and T, we average the temperatures corresponding to the same pixel 
in different matrix positions to have the observed temperature map Mt (including noises). 

In similar way, from the matrices N and W (or G) we obtain the observed temperature 
map M-w (or Mq) computed in presence of white noise only (or in absence of instrumental 
noises) . We have verified that the map Mg is identical to the "input" map for all the observed 
pixels, so confirming the validity of the geometrical part of our flight simulation code. 

5 Destriping techniques 

We have developed a technique in order to eliminate the effects of gain drifts in Planck 
signal due to the 1// noise effect. The method is derived from that proposed for the Phase 
A Study and re-analyzed by Delabrouille (1997). On the other hand, our treatment of the 
Planck observation simulation, although simplified, is general enough to be close to the "real" 
Planck observations; so we draw out here below the destriping mathematical formalism, in a 
way directly applicable to our simulated data streams. 

5.1 Mathematical formalism 

In this section we discuss how to eliminate the effects of gain drifts on timescales greater than 
that for which the spin axis points at a given direction (2 hours for the Phase A scheme, or 
1 hour or 40 minutes for possible new observational strategies in which the spin axis shift of 



2.5' or 1.6' in 1 hour or 40 minutes respectively), i.e. the satelhte scans a given circle in the 
sky. 

After we removed the drifts within any given scan circle by averaging the observations 
over the corresponding cycles (see section 3.1), each set of observations at a given spin axis 
direction, denoted by the index i, is characterized by an additive level Ai which is related to 
the "mean" 1// noise level during the observation in that scan circle. These levels Ai are 
different for different circles, due to gain fluctuations. Our goal is to obtain a reduced set of 
observations of different scan circles by removing the contamination that affects any circle. So 
we will subtract to all sets of observations on a given scan circle their own characteristic level 
Ai. As a variance, we can attribute more constant levels, say rii levels, per single scan circle. 
From the computational point of view this is exactly equivalent to rearrange all the matrices 
in section 3.3 by dividing their rows in ni parts that have to be appropriately relocated to 
construct new matrices with Ug x n/ rows and Up/rii columns that can be then analysed 
exactly as in the case of a single constant per circle. 

For estimating all these levels we use a computation scheme able to simultaneously find the 
pixels in common between different scan circles and generate a linear system whose solution 
gives the unknowns Ai. 

The observations of different directions in the sky explored by the satellite have been 
recorded in three matrices N, T, E of n^ rows and Up columns, where Ug is the number of 
different spin axis directions, and np is the number of samplings at different horn pointing 
directions in a given circle. Here Nn (z = 1, . . . , n^, / = 1, . . . , Up) contains the pixel number 
corresponding to the observed direction in the sky, and Tn and En are respectively the 
corresponding observed temperature (full signal due to the sky plus noises) and the estimates 
of the rms noise, essentially due to the white noise. We observe that En is properly related to 
the amplifier noise temperature and to the observed antenna temperature which depends in a 
real case on the local thermal conditions and on the true sky temperature but it may depend 
also on possible "spurious effects". We have in reality only "first-order" informations about 
all these quantities at this level; on the other hand its accurate knowledge is not crucial, being 
the amplifier noise temperature typically higher than the observed temperature. Anyway we 
hope to have more accurate informations from accurate thermal models and from iterating 
our data reduction scheme to have a good determination of sky temperature, which is of 
course our goal. 

In the following the first index will denote the row index and the second the column index. 

We must check for all the possible crossing points for two different circles for the whole 
ensemble of n^ circles (see also section 3.3). 

Let TT be an index that identifies a generic couple of different observations corresponding 
to the same pixel in the sky, i.e. a pixel in common between two scan circles: tt ranges from 
1 to ric, where Uc is the total number of couples found. Therefore the index vr is related to 
two elements in the matrix N: tt ^ {il,jm). Here i and j identify the two circles for which 
we found common pixels, / and m are the common pixels positions on the circle i and on 
the circle j respectively. So we have Nn = Njm- As a variance, we can replace the matrix 
N with one of the matrices Nh (see section 3.3 and 7.4) to search for the pixels in common, 
according to the adopted averaged maximum distance for recognizing two pixels in common. 

We want to minimize the quantity: 



all couples 



[{Ai - Aj) - {Til - Tj 



12 



jmj 



^il + ^"jm 



E 

7r=l 



[{A, - Aj) - {Til - Tjm)] 



(10) 



respect to the set of the unknown levels Ai] the index vr in the right hand side of this equation 
remembers that each set {il, jm) derive from a given pixel vr. S is quadratic in all the unknown 

10 



Ai] on the other hand, only the differences between the levels Ai enter in this expression, so 
that the solution will be indeterminate, i.e. the levels are determined apart from an arbitrary 
additional constant (with no physical meaning, as obvious for anisotropy measurements). 
To remove this indetermination, we add a constraint to the Ai quantities: 



h=l 



0. 



(11) 



This is equivalent to minimize the quantity: 



(12) 



\h=l 



Now let's go into some algebra. We perform the derivate of the previous equation, and 
finally we have: 



1 dS' 

2dA~k 



7r=l 



[{Ai - Aj) - {Til - Tjm)] ■ [Sik - Sjk] 



Efi + E' 



jm 



h=i 







(13) 



for all k = 1, . . . ,ns (here the 6 are the usual Kronecker symbols). So we have a set of Ug 
equations: 



J2CktAt = Bk, k=l,..., 



n^ 



(14) 



t=i 



We denote with C and B, respectively, the matrix of the coefficients Ckt and the vector of 
the coefficients B^. 

To be concrete, we show here as C and B are formed as we extract the pixels in common 
between the different rows. First of all, we set B = and C^t = 1 V/c, t (setting all C^t to 1 
takes into account the second term of eq. |l^). Then for each couple tt of pixels in common 
between two scan circle we define: 



X-n 



Efl + E'^jm 



and 



Til 



i-jm 



Efi + Ej^ 



(15) 



(16) 



Prom the above equation and the definition of Kronecker symbol, we easily have that a 
given couple vr contributes only to two equations of our linear system, those for k = i or 
k = j, where as usual i and j corresponds to two different observations of the same pixel. If 
we iteratively increment the coefficients of C and B as we find a new couple, explicitely we 
have [remember n — > {il,J7n)]: 



^ii *■ ^ii + XiT 



Gj' 



a 



v 



a 



J* 



c^ 



c^ 



33 

Bi 
Bi 



]i ^ Xtt 
Cjj + Xtt 
Bi +T^ 



B 



3 



(17) 
(18) 
(19) 
(20) 
(21) 
(22) 
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Summing up, we have that each couple vr contributes to only six terms, and the resulting 
system shows a complete symmetry with respect to the exchange of the indexes i and j. 

The linear system defined by eq. ( p!4[ ) has some interesting properties which considerably 
simplify the numerical computation of its solution. In particular, the matrix C: 

• is symmetric, so we can hold in memory only half of the matrix (say the upper-right 
part) and solve the system and speed-up the code by computing only half of the matrix 
coefficients. This is possible because the Gauss reduction algorithm preserves at each 
step the symmetry of remaining part of the matrix; 

• is positive defined, so we never find a null pivot when reducing a non-singular matrix 
(Strang 1976). This allow us to solve the system without having to exchange rows or 
columns, so preserving the symmetry; 

• is not singular (provided that there are enough intersections between scan circles), 
because the only fundamental indetermination has been removed by imposing the con- 
straint (|Tl|). 

Anyway, after the system has been solved, it is our care to replace the solution into the 
original C matrix, to verify its correctness and check for rounding errors and/or accidental 
degenerations. 

5.2 Some remarks on numerical efficiency, RAM requirements and off-sets 

In order to speed the construction of the matrix C and of the vector B, we have found that it 
is very advantageous to firstly order (we use the quick sort algorithm) all the elements of the 
matrix N, i.e. the observed pixels (4 bytes integers), in the first column of a new "matrix" U 
(of Up X rig rows and 3 "columns") by keeping memory of their locations (2x2 bytes integers), 
in the original matrix N in the other two "columns" of U. In this way we simply extract 
once for all each pixel in common between two scan circles for all scan circles, being the same 
pixel located in contiguous rows in the matrix U, by considering all the possible pairs of rows 
of U with the same element in the first column, with the simple caution that the elements 
of the second column of U, i.e. the original rows in the matrix N, are different. In this way 
the "scanning" of the matrix N and the construction of C and B according to the rules of 
section ^^ turns to be very fast. It is immediate to use the matrices Nh, containing the pixel 
numbers at higher resolutions, in the construction of C and B if one want to adopt more 
stringent conditions on the distance between pixels in common. In addition, working with 
the additional matrix U optimize the construction of the simulated maps from the simulated 
data streams (see section |^, being immediate to recognize in the matrix U when the same 
pixel has been observed. 

For the solution of the linear system ([l^ ) we have found that the Gauss elimination method 
works very well. We prefer to construct and solve the system by using double precision 
accuracy, to have high numerical accuracy and to be sure of avoiding artificial numerical 
singularity; also, due to matrix symmetry and positive definiteness, we do not need pivot. 

To build up the linear system, we have to keep the memory space for the system matrix, 
the system known terms, the observed temperature matrix and the auxiliary (integer) matrix 
U; by taking advantage of this symmetry, and by considering in general ni constant levels per 
scan circle the memory requirement is: Sbytes x [nins{nins + l)/2 + niUg + n^rip] + 4bytes x 
2 X UsUp. For example, at 30 GHz (FHWM~ 30', Up ~ 2100), for the case of 5' shift of the 
spin axis (n^ = 4320) we need about 220 (440) Mbytes by working with n/ = 1 (n; = 2); for 
a 2.5' shift of the spin axis, n^ = 8640 and we need about 590 (1500) Mbytes by working 
with ni = 1 {rii = 2). For sake of illustration, if we have a beam of ~ 10' (like the nominal 
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100 GHz beams) and we want to record 4 samplings per beam we will have Up ~ 8400; for a 
2.5' shift of the spin axis and by working with ni = 1 (ni = 2) the memory requirement is of 
about 1500 (2400) Mbytes. This memory problem can be solved by taking advantage of disk 
buffers; we discuss our solution in the Appendix B. 

For solving the system we only need to keep in memory the system matrix and known 
terms, and the memory requirement is: Sbytes x [nins{nins + l)/2 + niUs]- At 30 GHz and 
with 3 samplings per beam, we need about 75 Mbytes if Ug = 4320 and ni = 1, 300 Mbytes if 
Us = 4320 and n/ = 2 orn^ = 8640 and ni = I and 1200 Mbytes if n^ = 8640 and n/ = 2. This 
problem maybe crucial depending on the available amount of RAM, especially because the 
Gauss elimination continuously changes the system components. Also this memory problem 
can be solved by taking advantage of disk buffers (Appendix B). 

After the solution of the linear system, we obtain a "destripped" matrix D by subtracting 
the level Ai to the i— th row (i = l,ns) of the matrix T Then we apply to this matrix the 
same treatment of section Q and we have the observed destripped temperature map Md- 

We observe that, contrary to the case of pure white noise, the average of a l//-type noise 
series can be significantly different from zero. Then, the map Mt (as well as the map Md, 
but in general with a somewhat different value) may present an off-set with respect to the 
map Mg; on the contrary the off-set between the map Mw and the map Mq is negligible. 
As a typical example, for the simulation with a = 90° (see Table 1) we find that these off-sets 
are ~ 4.8^K; for comparison, the off-set we find between the maps Mw and Mg is much 
smaller, ~ 0.075^K. The off-set between Mt (Md, Mw) and Mq is of course not relevant 
for anisotropy measurements, nor we are able to subtract it in a real case (on the other hand 
we must pay attention to the fact that off-sets may be present between maps produced by 
different receivers at the same frequency). The off-set of the map Mt (Md, Mw) must be 
removed by subtracting the difference between the average of the map Mt (Md, Mw) and 
of the map Mq: this is necessary for a correct quantitative analysis of stripes magnitude and 
destriping efficiency. 

We indicate with Mt, Md (and Mwi but it is not relevant in practice for this matrix) 
the above maps, when this kind of off-set has been removed. 

6 Estimators of the destriping efficiency 

In the previous sections we have described our simulations of the Planck observations and 
the basic treatment to convert observational data streams into sky maps, including destriping 
techniques. Here we analyse the efficiency of the adopted destriping technique, by considering 
well known estimators. 

6.1 Ratio between the x^'s 

We expect that the average of the squares of differences between the elements of the maps 
Mw and Mg divided by the observation sensitivity [essentially the estimator Xr'j '^^ ^i^^ 
call it in this case (Xr)w^] is very close to 1, because only the white noise is present in this 
case. We can compute the same estimator for Mt — Mg [{Xr)T, undestripped case] and 
Md — Mg [{Xr)D destripped case]. (For the above consideration, by using observed maps 
without removing the off-sets one will find a meaningless amplification of the Xr)- 

On the other hand, we find that the exact value of {Xr)w may be just a little different 
from 1 depending indeed on the assumed sensitivity; for example it is just a little different is 
we divide Mw — Mg by the map Mg (i.e. by using the sensitivity proper of any pixel - we 
will use this definition in our tables) or the average of the sensitivities in the map Mg or the 
estimate of the average sensitivity obtained on the basis the global mission time, the observed 
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number of pixels (393216 for our maps at COBE-cube resolution 9) and the properties of 
considered receiver. Then we prefer to use the ratio between the {Xr)T ^ or {Xr)D ~ and 
{Xr)w as estimator. This "renormalized" estimator (that we will denote i>y xf-nT ^^'^ XrnD 
respectively for the undestripped and destripped cases) results independent of the choice of the 
sensitivity adopted for the Xr calculation, so allowing a better understanding of the magnitude 
of striping effect and of destriping efficiency. We will quantify the destriping efficiency by using 
the relative decrease of the renormalized Xr^ i-6- with the quantity [(Xr)-D~(Xr)T]/[(Xr)T — !]• 

6.2 Magnitude of stripes temperature 

From the values of XrnT ^^'^ Xrn D defined above and from the average rms white noise, 
rmsw, for the observed pixels derived from the sensitivity map Mg we can easily give an 
estimate of the rms temperature of the stripes before, rmsT, and after destriping, rmso- 
Under the hypothesis that the error introduced by the noises in each pixel may be thought 
as a sum of two uncorrelated contributions from white noise and 1// noise, we have rmsT = 

rmswJxr,n,T ~ 1 ^'^^ rmsD = rmswJxr,n,D ~ 1 [Method (a)]. 

Another estimate [Method (b)] of the rms stripes temperature can be obtained by directly 
evaluating the global temperature rms difference before, rmstot,T, or after, rmstot,D, the 
destriping from the comparison with the maps Mq and by assuming that they are given by 
the sum in quadrature of rmsw and rmsT or rmso. From the values of rmsw, T~mstot,T and 
rmstot,D we can calculate rmsT or rmso- 

Finally [Method (c)], we can treat the 1// contribution to the total noise like a systematic 
(and not statistical) error and therefore to assume that rmstot,T or rmstot,D are simply given 
by the sum of rmsw and rmsx or rmsD- 

For estimating the destriping efficiency in terms of residual stripes temperature we will by 
using the relative decrease of the stripes temperature, i.e. the quantity {rmsD — rmsT)/rmsT- 

7 Results 

From the visual inspection of the simulated maps the effect of the noises is of course not 
clearly evident; we only recognize a somewhat degradation of the map details. The stripes 
become more evident when we plot the noises map only (see Figure 3); they can be obtained 
by subtracting the CMB fiuctuation map (Mq) to the observed map. The stripes figures 
well reproduce the adopted scanning strategy. 

7.1 Destriping versus scanning strategy 

The visual inspection of our "stripes" maps do not allow to quantify the striping magnitude 
and its reduction obtained from the destriping procedure (see Figures 3 and 4). 

The statistical analysis of the maps allows a much better understanding of the destriping 
procedure efficiency. We present here (see Tables 1-^4) the results of our simulations (for a 
typical channel at 30 GHz) in terms of reduced Xr n and of stripes temperature for the three 
choosen values of the (constant) angle a between the spin axis and the telescope direction 
and for the considered case with a = 90° and spin axis oscillations (see also the Appendix 
A). We report also in the tables some informations on relevant quantities: the ideal white 
noise for a single sampling time which weakly increases with a for geometrical reasons if we 
want to have 3 samplings per beam; the (single receiver) average rms white noise, rmsw, 
(expressed in mK) for the observed pixels for a 360 days mission; the square of the ratio R 
between these white noises, normalized to the intermediate case of a = 85°; the percentage of 
sky which results to be observed by the considered single off-axis beam, which also increases 
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with a; the number of couples in common find in the destriping procedure and the number 
of constant levels per scan circle used in the destriping code; the map resolution and the 
resolution used for searching the pixels in common. We remember the adopted value, 0.05 
Hz, of 1// knee frequency and of the bandwidth, 6 GHz, and that the sampling time is of 
about 0.03 sec, the exact value depending on the choosen value of a. In the tables we report 
our values of Xrn before and after the destriping procedure and the relative (%) decrease 
without the multiplicative factor R^ and by taking it into account. Indeed the (white noise) 
sensitivity per pixel is different for different scanning strategies; then, by including the factor 
R^ we "renormalize" the values of Xr n ^^ ^^^ same (white noise) sensitivity level, so making 
the them essentially independent of the a-dependent sensitivity. 



Table 1: destriping results; a = 90°. 



Some global parameters 


ATw = 1.322mK 
Map Res. = 9 


rmsw = 27.18^K 
Res. common pix. = 9 


i?2 = 1.00666 
ni = l 


pix. 


% sky = 99.98 
in common = 2.70 x 10^ 


Before destriping 


After destriping 


% improvement 




Method 


xln,T = 1-1770 
xln,TR'' = 1-1848 
rmsT = 11.4^K 
rmsT = 17.3;uK 
rmsT = 5.03^K 


xl^^D = 1-0216 
xln,DR' = 1-0284 
rras£, = 3.99^K 
rrasD = 6.93^K 
rmsD = 0.87^K 


87.8 % 

84.6 % 
65.1 % 

59.9 % 

82.7 % 




(a) 
(b) 
(c) 



Table 2: destriping results; a = 85°. 



Some global parameters 



ATw = 1.318mK 
Map Res. = 9 



rmsw = 27.09^K 
Res. common pix. = 9 



i?2 



% sky = 99.43 
pix. in common = 2.34 x 10^ 



Before destriping 



After destriping % improvement 



Method 



A/2 
Xr,n,T - 


-- 1.2709 


xln,D = 1-0142 


94.7 % 


xin.T^' 


= 1.2709 


xln,DR' = 1-0142 


94.7 % 


rmsT = 


14.1^K 


rmsD = 3.23^K 


77.1 % 


rmsT = 


21.6^K 


rmsD = 7.31^K 


66.2 % 


rmsT = 


7.98^K 


rmsD = 0.969^K 


87.9 % 



(a) 
(b) 
(c) 
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Table 3: destriping results; a = 80°. 



Some global parameters 


ATvK = l.SllmK 
Map Res. = 9 


rmsw = 26.90^K 
Res. common pix. = 9 


i?2 = 0.98602 
ni = 1 


pix. 


% sky = 98.12 
in common = 2.20 x 10^ 


Before destriping 


After destriping 


% improvement 




Method 


xln^T = 1-0396 
xln,TR' = 1-0251 
rmsT = 5.35^K 
rmsT = 8.05fiK 
rmsT = 1.1 7^K 


xln,D = 1-0292 
X?,„,D«' = 1-0148 
rmsD = 4.60^K 
rrasD = 6.94^K 
rmsD = 0.88^K 


26.19 % 

41.0 % 

14.1 % 
13.8 % 
24.8 % 




(a) 
(b) 
(c) 



Table 4: destriping results; a = 90° it 10° (10 sinusoidal oscillations). 



Some global parameters 



ATw = 1.322mK 
Map Res. = 9 



rmsw = 29.07^K 
Res. common pix. = 9 



R^ = 1.15152 % sky = 100 

ni = 1 pix. in common = 2.36 x 10^ 



Before destriping 



After destriping 



% improvement 



Method 



A/2 
Xr,n,T - 


-- 1.1165 


v2 


-- 1.0153 


86.9 % 


xin.Tli' 


= 1.2857 


Xr,n,D^'^ 


= 1.1691 


40.8 % 


rmsT = 


9.92^K 


rmsD = 


3.60^K 


63.8 % 


rmsT = 


15.2^K 


rmsD = 


8.62^K 


43.3 % 


rmsT = 


3.73^K 


rras£, = 


1.25^K 


66.5 % 



(a) 
(b) 
(c) 



7.2 Destriping versus 1// knee frequency 

For the interesting case a = 90° with no oscillations, we carried out other a simulation with 
a much larger value of the 1// knee frequency, /^ = lOHz, of order of that expected for total 
power radiometers. 

It is interesting to study the stripes effect and destriping performance under this very 
pessimistic condition. 

Tables 5 shows our results that have to be compared with those of Table 1, based on the 
theoretical estimate of fk of our kind of radiometers. 
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Table 5: destriping results; a = 90°; fk = lOHz. 



Some global parameters 


ATvK = 1.322mK 
Map Res. = 9 


rmsyi/ = 27.18^K 
Res. common pix. = 9 


R^ = 1.00666 

ni = 1 


pix. 


% sky = 99.98 
in common = 2.70 x 10^ 


Before destriping 


After destriping 


% improvement 




Method 


xln,T = 9-2846 
X?,n,T^' = 9-3464 
rmsT = 78.2^K 
rmsT = 258^K 
rmsT = 233^K 


xln,D = 2-6769 
xIu^dR^ = 2.6947 
rmsD = 35.4^K 
rrasD = 67.7^K 
rmsD = 45.8^K 


79.8 % 

79.7 % 

54.8 % 
73.7 % 
80.3 % 




(a) 
(b) 
(c) 



7.3 Destriping with more than one constant per scan circle 

For the reference case a = 90°, both with the theoretical prediction for /^. and for the case 
with fk representative of total power radiometers, we have applied our destriping code by 
using two constants per scan circle. The results are shown in Tables 6 and 7 that must be 
compared with Tables 1 and 5 respectively. We find that the use of more constant per circle 
does not help the destriping technique. 



Table 6: destriping results; a = 90° 



Some global parameters 


ATw = 1.322mK 
Map Res. = 9 


rmsw = 27.18^K 
Res. common pix. = 9 


R^ = 1.00666 

ni = 2 


pix. 


% sky = 99.98 
in common = 2.70 x 10^ 


Before destriping 


After destriping 


% improvement 




Method 


xln,T = 1-1770 
xI^.^tR'' = 1-1848 
rmsT = 11.4^K 
rmsT = 17.3^K 
rmsT = 5.03^K 


xln,D = 1-0280 
xI^^dR'' = 1-0348 
rmsD = 4.55//K 
rmsD = 7.62^K 
rmsD = l.OS^K 


84.2 % 
81.2 % 

60.1 % 
56.0 % 

79.2 % 




(a) 
(b) 
(c) 
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Table 7: destriping results; a = 90°; fk = lOHz. 



Some global parameters 


ATvK = 1.322mK 
Map Res. = 9 


rmsyi/ = 27.18^K 
Res. common pix. = 9 


R^ = 1.00666 

ni = 2 


pix. 


% sky = 99.98 
in common = 2.70 x 10^ 


Before destriping 


After destriping 


% improvement 




Method 


xln,T = 9.2846 
X?,n,T^' = 9.3464 
rmsT = 78.2^K 
rmsT = 258^K 
rmsT = 233^K 


xln,D = 2.7584 
xIu^dR^ = 2.7768 
rmsD = 36.2^K 
rrasD = 70.0^K 
rmsD = 47.9^K 


78.8 % 
78.7 % 
53.7 % 

72.9 % 
79.4 % 




(a) 
(b) 
(c) 



7.4 Destriping versus distance conditions 

For the reference case a = 90°, both with the theoretical prediction for /^ and for the case 
with /fc representative of total power radiometers, we have applied our destriping code by 
using the map pixels at higher resolution to search for pixels in common. The results are 
shown in Tables 8 and 9 that must be compared with Tables 1 and 5 respectively. We conclude 
that the use of a more stringent condition to find the coincidences of the pointing directions 
in different scan circles does not help the destriping technique. 



Table 8: destriping results; a = 90°. 



Some global parameters 



ATw = 1.322mK 
Map Res. = 9 



rmsw = 27.18;^K 
Res. common pix. = 10 



R^ = 1.00666 

ni = 1 



% sky = 99.98 

6.97 X lO'^ 



pix. m common 



Before destriping 



After destriping 



Vo improvement 



Method 



A/2 
Ar,n,T - 


-- 1.1770 


xl,n,D = 1-0235 


86.7 % 


A/2 p2 


= 1.1848 


xln,DR' = 1-0303 


83.6 % 


rmsT = 


11.4^K 


rmsD = 4.73/^K 


58.5 % 


rmsT = 


17.3^K 


rras£, = 7.16/uK 


58.6 % 


rmsT = 


5.03^K 


rmsD = 0.928;[iK 


81.6 % 



(a) 
(b) 
(c) 
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Table 8': destriping results; a = 90° 



Some global parameters 


ATvK = 1.322mK 
Map Res. = 9 


rmsw = 27.18/uK 
Res. common pix. =11 


R^ = 1.00666 

ni = 1 


pix. 


% sky = 99.98 
in common = 1.62 x 10^ 


Before destriping 


After destriping 


% improvement 




Method 


xln,T = 1-1770 
xln,TR' = 1-1848 
rmsT = 11.4^K 
rmsT = 17.3;uK 
rmsT = 5.03^K 


Xl^^n = 1-0334 
xln,DR' = 1-0403 
rmsD = 5.46/LiK 
rmsD = 8.10/[/K 
rmsD = 1.18/uK 


81.1 % 

78.2 % 

52.1 % 

53.2 % 
76.5 % 




(a) 
(b) 
(c) 



Table 9: destriping results; a = 90°; fk = lOHz. 



Some global parameters 


ATw = 1.322mK 
Map Res. = 9 


rmsw = 27.18/uK 
Res. common pix. = 


10 


R^ = 1.00666 
ni = 1 


pix. 


% sky = 99.98 
in common = 6.97 x 10^ 


Before destriping 


After destriping 




% improvement 




Method 


xln,T = 9-2846 
Xr,n,T^' = 9.3464 
rmsT = 78.2/iK 
rmsT = 258^K 
rmsT = 233^K 


xln^D = 2-6970 
xIu^dR^ = 2.7150 
rmsD = 23.0^K 
rmsD = 68.3/nK 
rras£, = 46.3/uK 




79.5 % 

79.5 % 

70.6 % 
73.5 % 
80.1 % 




(a) 
(b) 
(c) 
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Table 9': destriping results; a = 90°; fk = lOHz. 



Some global parameters 


ATvK = 1.322mK 
Map Res. = 9 


rmsyi/ = 27.18/uK 
Res. common pix. =11 


R^ = 1.00666 

ni = 1 


pix. 


% sky = 99.98 
in common = 1.62 x 10^ 


Before destriping 


After destriping 


% improvement 




Method 


xln,T = 9-2846 
xln^rR^ = 9-3464 
rmsT = 78.2^K 
rmsT = 258^K 
rmsT = 233^K 


xln,D = 2-7752 
xIu^dR^ = 2-7937 
rmsD = 24.2^K 
rmsD = 70.5/uK 
rmsD = 48.4;uK 


78.6 % 
78.5 % 
69.0 % 

72.7 % 
79.2 % 




(a) 
(b) 
(c) 



8 Discussion and conclusions 



An analytical estimate of the maximum excess noise factor, F, due to the stripes related to 1// 
effect has been given by Janssen et al. (1996); they found F ~ [l + r/fe(21nnp + 0.743)]-*^/^. By 
thinking F^ as equal to l + {l^rms/rmsw)'^, we have a fractional additional rms noise respect 
to the rms noise, rmsw-, obtained in the case of pure white noise given by {d^rms / rmsy/)'^ = 
Tfk{2\nnp + 0.743). For example for our simulations we have Up ~ 2100, a sampling time 
r ~ 28msec, corresponding to an angle of ~ 10' in the sky, and fk = 0.05Hz (or 10 Hz) . 
With these number we get {Arms/rmsw)^ — 0.022 (or 4.5) for a pixel of 10' x 10'. Our map 
pixel is 19.4' x 19.4' (COBE-cube resolution 9); therefore rms^^ reduces by a factor ~ 4 and 
we expect to have an additional (Arms / rmsw)'^ per pixel, i.e. an additional reduced x^i 
of about 0.1 (or 18). The results shown in our tables are in quite good agreement (always 
within a factor 2) with these analytical estimates. Somewhat larger values may be expected 
from the larger observational time toward high ecliptic latitudes, the consequent reduction 
of the white noise and the increasing of the relative weight of the 1// noise. We stress here 
that, contrary to the case of pure white noise, even for the same scanning strategy the final 
effect of 1// noise may be quite different for different simulations: larger or smaller effects 
can be obtained according to the (simulated) behaviour of 1// gain fluctuations. Only from 
a very large set of simulations for each considered scanning strategy we can derive a robust 
evaluation of the "averaged" final effect of 1// noise (see below). For this reason, the results 
shown in our tables have to be considered as first order estimates of the final 1// noise effect 
expected in a given scanning strategy rather than detailed predictions. 

Without applying the destriping procedure, the typical amount of stripes temperatures 
(i.e. the rms values) per resolution element (of ~ 19.4' x 19.4') ranges from few ^K to one 
or two tens of fiK (according to the method adopted for estimating them) and is always 
significantly less than the corresponding single beam sensitivity of about 27 -^ 30/iK. 

The residual stripes that remain after the reduction through our destriping code show 
typical rms temperatures of few //K, roughly independently of the adopted scanning strategy. 
The efficiency of the destriping algorithm is quite good ('-^ 20% -^ 90%, according to the 
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adopted estimator and depending in part on the scanning strategy). We stress that in any 
case, for our kind of radiometers, the reduced additional noise is of few percent in terms of 
increased reduced x^ £^^d of few //K in terms of stripes temperature. This is not particularly 
critical. Nevertheless it must be compared with the final Planck sensitivity as it results by 
combining the data from all the receivers at the same frequency (see below). 

In the case of much larger 1// contaminations like those expected for higher values of 
the knee frequency in the case of total power radiometers, the final impact on observed sky 
maps is much higher. The destriping code allows to reduce a large fraction of the added noise 
(see Table 5), nevertheless a significant increase of the final noise still remains. Then, the 
importance of the reduction via "hardware" of the 1// noise is strongly recommended. 

About the dependence of the striping effect and of the destriping efficiency on the scanning 
strategy, our preliminary results suggest the following conclusions. 

• The magnitude of the sensitivity degradation due to the stripes is only weakly dependent 
on Q for undestripped maps as well as for destripped maps (see Tables 1-^3). 

• The final residual added noise is roughly independent of the scanning strategy and on 
the unreduced added noise, almost not relevant unreduced contaminations (see Tables 
1-4). 

• Oscillations of the spin axis do not improve significantly the destriping efficiency (see 
Table 4). On the contrary it is well known that oscillations of the spin axis may 
introduce further systematic effects related to variations of the illumination by the sun 
and of shielding performances. 

• Our results are of course related to the beam position in the sky field of view. For 
on-axis beams and for beams located at (pB ~ 0° the case at a = 90° is expected to 
give the worst results, whereas even in this case the destriping results are expected to 
improve for off-axis beams located at (/)b ~ 90°. 

• We find that, in spite of the larger computing time, the efficiency of the destriping 
technique does not improve by using two level constants per scan circle (see Tables 6 
and 7) . Indeed the number of "physical" conditions (the number of pixels in common) 
do not depend on the choosen number of level constants adopted in the destriping 
procedure. Then, by searching (in the solution of the linear system, see section 5.1, 
obtained by the condition of minimization of S) for a number of unknown larger by a 
factor two, by using the same number of informations, we expect that the uncertainty 
of each unknown will be larger. We infer that the advantage of using more constant per 
scan circle found by Delabrouille (1997) for a similar knee frequency but for the case of 
thermal drifts (noise spectrum proportional to 1//^) has to be related to the different 
kind of noise spectrum. The noise fluctuations on long timescales will be higher in the 
case of 1//^ noise spectrum than in the case of 1// noise spectrum; in the former case 
the use of more constants per circle may allow a more appropriate subtraction of the 
gain fluctuations in each circle, whereas in our case this advantage is balanced by the 
increased uncertainty in the determination of the levels and the global effect results in 
to a small decreasing of destriping efficiency. 

• The use of more stringent conditions for identifying the pixels in common does not 
improve (see Tables 8 and 9) the destriping results. This option allows to reach a more 
accurate superposition of the pixels used in the destriping; nevertheless the number of 
couples decreases about by a factor 4 (or 16) by using pixels 4 (or 16) times smaller, 
and so a number of conditions significantly smaller than that found before enters in the 
minimization of S (see section 5.1). 
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Given all these results we believe that a simple scanning strategy with constant a in the 
range 85° -^ 90° can offer the advantages of a large (practically full) sky coverage, quite good 
destriping performances together with the minimization of all the systematic effects related 
to the variations of thermal conditions. 

We draw here below a brief guide-line for the future simulation work on this topic. 



• 



We intend to test different sampling strategies: for example with a smaller shift (2.5') of 
the spin axis direction and a corresponding smaller observation time for spin axis direc- 
tion (1 hour). Such a kind of sampling, possibly with the same (~ 8700, corresponding 
to 4 samplings per beam at 100 GHz - FWHM~ 10') number of samplings per scan 
circle independently of the frequency, may be the final LFI/Planck sampling. A strat- 
egy of this kind allows to have a larger number of pixels in common in the destriping 
procedure; on the other hand the sensitivity of the temperature measurements for the 
different scan circles degrades. It may be interesting to investigate pro and contro of a 
variance of this kind from the point of view of the 1// noise reduction. 

We may be also interested in an accurate verification of the validity of the choosen 
telescope rotation velocity around the spin axis. 

It may be interesting to apply our codes to higher ratios between the beam FWHM 
and the map pixel size, for example by working, at 30 GHz, with maps at COBE- 
cube resolution 10. Indeed, considering input maps at a higher resolution roughly tests 
the importance of assuming a better efficiency in the data streams deconvolution, for 
example by fully exploiting the beam oversampling. 

As better approximation, it is interesting to implement the convolution with the beam 
for a moving telescope in the observation simulation code and to search for robust 
and fast criteria for establishing in this situation when it is possible to consider that 
integrations in different scan circles can be really referred to the same direction in the 
sky. Of course this problem is correlated to the technique adopted for deconvolving the 
data streams to obtain observed maps; on the other hand, the stripes magnitude must 
be small enough to not significantly alter the deconvolution procedure. 

We have found that the noise added by the 1// effect is not particularly critical, com- 
pared to the single beam (white noise) sensitivity. Nevertheless it must be compared 
with the final Planck sensitivity as it results by combining the data from all the re- 
ceivers at the same frequency. Indeed, we do not expect that the 1// noise magnitude 
decreases as the square root of the number of receivers, as white noise does: by carrying 
out several simulation with the same set of physical parameters for the same scanning 
strategy and averaging the corresponding maps, we can address this topic. 

We intend to apply the methods of inversion of CMB maps (Muciaccia et al. 1997) for 
deriving the angular power spectrum of the observed maps. By comparing the CMB 
angular power spectrum obtained in presence of 1// noise contamination with that 
derived in the case of pure white noise (and of course with that of the input map) it 
is possible to estimate the 1// noise impact on the extraction on the key cosmological 
informations, almost in the case of gaussian fluctuations like those expected in infla- 
tionary scenarios. Particular attention has to be attempted for evaluating the impact 
on our science in the context of topological defects, like cosmic strings for example, 
which introduce non gaussian features in the power spectrum. The characteristic geo- 
metrical pattern of 1// noise stripes (related to the scanning strategy) have to be used 
for disantangle between instrumental and cosmological deviations from the gaussianity. 
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• Of course, we intend to extend in the next future our analysis to Planck measurements 
at higher frequencies. 

The full success of missions like Planck and MAP require a good control of all the relevant 
sources of systematic effects. Discrete sources above the detection limit must be carefully 
removed and accurate models for foregrounds radiation and anisotropies (Brandt et al. 1994, 
Danese et al. 1996, Bouchet et al. 1997, Toffolatti et al. 1995, 1997) are required to keep 
the sensitivity degradation in the knowledge of CMB anisotropies below few tens percent 
(Dodelson 1997). Optical distortions, which produce a non-symmetric beam response for 
feed-horns located away from the centre of the focal plane, introduce other systematic effects; 
they must be minimized by optimizing the telescope and the focal plane assembly design 
(Mandolesi et al. 1997b). Thermal drifts (Bersanelli et al. 1996), which couple to the 1//- 
type noise here discussed, can also generate stripes in the observed maps; efficient shield is 
required together with accurate reduction of sidelobe effects and optimization of the thermal 
conditions during the mission. 

All in all maximum efforts should be addressed to optimize the cooling efficiency and the 
observational strategy and to improve the methods for the data analysis in order to reduce 
the magnitude of the striping effect and of the other instrumental systematic effects. 
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Appendix A: Geometrical transformations between coordinate 
systems 

Let i,j,k the standard unit vectors in ecliptic coordinates and s a unit vector along the 
satellite spin axis outward the Sun direction. Let i the angle between s the and k (i.e. the 
ecliptic colatitude of s) and (p the angle between i and s (i.e. the ecliptic longitude of s) . For 
general scanning strategies i will be described by a, possibly not constant, function i = i{(p). 
Then s = sin icos (p i + sin isin (j) j + cos i k. Let i' = s and k' a unit vector hortogonal to s on 
the plane identified by the vectors k and s, namely k' = — cosicos^ i — cosisinc/) j + sini k. 
Let j' = k' Ai' (here A indicates the vector product). Let p the unit vector that identifies 
the pointing direction of the telescope optical axis. In the reference i' ,j' , k' the vector p can 
be defined by two angles: the angle a from s {a = 70° for the Phase A study, Bersanelli 
et al. 1996) and the angle, ip, between its projection on the plane identified by j', k' and 
k', with the convention p= cos a i' + sinasin^ j' + sinQcosV' k' . Given i' ,j' , k' in terms of 
i-ij^k it is easily to derive p in the same basis. We choose two coordinates x and y on the 
plane tangent to the celestial sphere in the telescope optical axis direction, p, with unit vector 
u and V respectively; we choose the x axis according to the condition that the unit vector 
u points always toward the satellite spin axis; indeed, for standard Planck observational 
strategy, this condition is preserved as the telescope scans different sky regions. With this 
choice of reference frame, we have that v = p f\ s/\p f\ s\ and u = v f\p/\v f\p\. In general, the 
coordinates (a:o,yo) of the beam centre in a ("satellite") reference XT,yT,ZT, corresponding 
to the unit vectors n, v,p, can be identified by two angles; we use here the colatitude 9b and 
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the longitude (j)B in this reference. Finally, the pointing direction of this generic (on-axis or 
off-axis) beam is given by the unit vector B = cos 9b p + cosc/i^sin^^ u + smcpBsinOB v. 

For sake of illustration, we consider here briefly three different kinds of scanning strategies 
(the first two options have been considered in the present simulations). 

• Spin axis always on the ecliptic plane. 

In this simple case we have i = 90° and s = cos (j) i + siiKJ) j. 

Given the angle a, we need to give the time dependences, cp = <j){t) and ip = tpit), of 
the spin axis longitude and of the telescope projection to fully determine the scanning 
strategy. The reference case is to change (j) of a certain angle Acj) (5', for example) after 
a given time interval (2 hours, for example) and to choose a given spin frequency (1 
r.p.m., for example) of the continuous rotation of ■0- 

• Sinusoidal oscillations of the spin axis. 

In this case we need also to define the amplitude of the oscillations, S (10°, for example), 
and the number of complete oscillations, nose (10 oscillations for 360 days, for example), 
per a complete rotation of the spin axis over the ecliptic (in 360 days, for example). 
Then, by choosing i = when i^ = 0, we have simply i = 6sm{nosc4') ■ 

• Precession of the spin axis. 

This case is just a little more complicate. We can consider a further unit vector / 
which moves always on the ecliptic plane; its ecliptic longitude, r] = r]{t) defines it: 
/ = COST/ i + sin?7 j. Let j" = — sinry i + cosry j. The satellite spin axis s precedes 
around / (we can consider for example again 10 precessions for 360 days, per a complete 
rotation of the axis / over the ecliptic plane in 360 days). Let ^ the angle between its 
projection on the plane identified by j" and k and the vector k and 6 the angle between 
/ and s (10°, for example). Then, the relation s = cos 5 f + sinJcos.^ k + sinJsin^ j" 
easily gives the colatitude i and the longitude (p of spin axis. 

Appendix B: System creation and solution with low memory 
usage 



As discussed in section |5.2| , the creation of the linear system requires a large amount of 
memory, usually much more than the available RAM. This problem can be avoided by taking 
advantage of disk buffers, essentially by splitting a large matrix into smaller blocks and 
creating it a block at a time. Our strategy to do this is very simple: 

1. a memory buffer is created, large enough to keep L lines; 



2. the algorithm described in section 5.1 is performed: for each couple vr of pixel in common 



between two scan circles i and j, the quantities Xn fmd Tj^ are evaluated; 

3. if i is in the range [0, . . . , L — 1], then equations (17), (18) and (21) are applied; if j is 
in the range [0, . . . ,L — 1], then equations (19), (20) and (22) are applied; 

4. after all the couples have been evaluated, the memory buffer is saved in a file; steps 
2 — 4 are then repeated for i and j in the range [L, . . . , 2L — 1], and then for i and j in 
[2L, ... , 3L — 1], and so on until all the A^ matrix lines have been created. 
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At the end of this loop, the Unear system has been created and stored in a file in binary 
format; the coefficients are organized in a matrix of A^ rows and A^ + 1 columns, where the 
last column contains the known terms. 

It is easy to see that this strategy saves memory space but increases considerably CPU 
time, because time required to create the whole matrix is directly proportional to the number 
of pieces N/L into which the linear system is divided. Due to the fact that program execution 
time is essentially dominated by the routines which solve the system, while creation time is 
negligible for the present practical purpose, we didn't bother to improve our code. 

A very different situation occurs for system solution, because both a large amount of 
memory and a lot of computation time are usually required. While time efficiency can't be 
usually improved (except in a few particular cases, when the matrix which defines the system 
has special symmetry properties and a lot of zero coefficients), memory requirements can 
be considerably reduced by taking advantages of disk swapping. However, a special care is 
required when choosing the strategy for disk operations, because disk time can easily blow 
up and overtake CPU time. Our algorithm behaves as follows: 

1. as for system creation, a memory buffer is allocated, large enough to keep L lines; 

2. then the first L lines (0, . . . , L — 1) are loaded in memory, and complete Gauss elimi- 
nation is performed on them: each line is reduced by the preceding lines and is used to 
reduce the following lines; 

3. each of the remaining N — L lines is sequentially load into memory, and reduced by 
each of the L lines stored in the memory buffer. In this way, we perform L steps of the 
Gauss elimination algorithm with a single disk operation; 

4. the memory buffer is flushed, and the next L lines {L,...,2L — 1) are loaded into 
memory, reduced by one another and used to sequentially reduce all the remaining 
lines of the system. Steps 3-4 are repeated until the original matrix has been completely 
reduced. 

As can be seen, our algorithm differs from the "standard" version only in the order in 
which the elimination is performed; in this way, we can limit memory requirements without 
increasing total solution time. It is easy to see that the total time spent for disk operation is: 

idisk = /3y ( - + 1 j (AT + 1) « /3disk-^ (23) 

CPU time scales as the cube of the linear size of the system: 

tcpu = ^cpuN^ (24) 

So we conclude that the ratio rj between disk time and CPU time is independent from 
the size A^ of the linear system, and only depend on the total number of lines L then we can 
hold simultaneously in memory: 

tdisk Pdisk /oir\ 

^ = -— = -^ — :j (25) 

Of course, the memory necessary to hold a line is proportional to the line length, so we 
need a larger amount of memory to solve a larger system; however, the size of the buffer is 
linear in system size A^, while the memory required to hold all the coefficients simultaneously 
is quadratic in A^. 
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FIGURE CAPTIONS 



Figure 1: The simulated (input) map of CMB anisotropies (CDM model, the dipole term 
is neglected). Galactic coordinates have been used for the plot. 

Figure 2: Schematic representation of the observational geometry. 

Figure 3: The unreduced simulated noise map (white plus 1// noise) for the simulation 
with a = 85°. Note the two small circular regions close to the ecliptic poles that are not 
observed by the considered off-axis beam; for graphic purposes we have filled them with a 
random noise distribution with variance given by the noise variance of the observed pixels. 
Note also the elongated sky region with noise significantly larger than the average, which 
corresponds to the sky regions that are observed a single time only, due to the choosen value 
of a. (Galactic coordinates have been used for the plot). 

Figure 4: The reduced simulated noise map (white plus "reduced" 1// noise) for the sim- 
ulation with a = 85°. Note again the two small circular regions close to the ecliptic poles 
filled with a random noise distribution. We note that the the elongated sky region with noise 
significantly larger than the average disappears as result of the destriping procedure. Also 
the stripes in the sky became much less evident. (Galactic coordinates have been used for 
the plot). 
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